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ABSTRACT 


We present the analysis of seven microlensing planetary events with planet/host mass ratios q < 10-4: 
KMT-2017-BLG-1194, KMT-2017-BLG-0428, KMT-2019-BLG-1806, KMT-2017-BLG-1003, KMT-2019- 
BLG-1367, OGLE-2017-BLG-1806, and KMT-2016-BLG-1105. They were identified by applying the Ko- 
rea Microlensing Telescope Network (KMTNet) AnomalyFinder algorithm to 2016-2019 KMTNet events. A 
Bayesian analysis indicates that all the lens systems consist of a cold super-Earth orbiting an M or K dwarf. To- 
gether with 17 previously published and three that will be published elsewhere, AnomalyFinder has found a total 
of 27 planets that have solutions with q < 1074 from 2016-2019 KMTNet events, which lays the foundation 
for the first statistical analysis of the planetary mass-ratio function based on KMTNet data. By reviewing the 27 
planets, we find that the missing planetary caustics problem in the KMTNet planetary sample has been solved 
by AnomalyFinder. We also find a desert of high-magnification planetary signals (A = 65), and a follow-up 
project for KMTNet high-magnification events could detect at least two more q < 1074 planets per year and 
form an independent statistical sample. 
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1. INTRODUCTION Among current exoplanet detection methods, a unique ca- 

pability of the gravitational microlensing technique (Mao 

& Paczynski 1991; Gould & Loeb 1992) is to detect low- 

Corresponding author: Weicheng Zang mass (Myianec S 20M) cold planets beyond the snow line 
3130102785 Q zju.edu.en (Hayashi 1981; Min et al. 2011), including Neptune-mass 
cold planets, which are common (Uranus and Neptune) in 
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our Solar System and cold terrestrial planets, which are ab- 
sent in our Solar System. Because the typical host stars of 
the microlensing planetary systems are M and K dwarfs, de- 
tections of q < 1074 planets (where q is the planet/host mass 
ratio) can reveal the abundance of low-mass cold planets and 
answer how common the outer solar system is. 

However, since the first microlensing planet, which was 
detected in 2003 (Bond et al. 2004), the first 13 years of mi- 
crolensing planetary detections only discovered six q < 1074 
planets! and none of them had mass ratios below 4.4 x 1075. 
The paucity of detected q < 1074 planets led to impor- 
tant statistical implications for cold planets. Suzuki et al. 
(2016) analyzed 1474 microlensing events discovered by the 
Microlensing Observations in Astrophysics (MOA) survey 
(Sako et al. 2008) and formed a homogeneously selected 
sample including 22 planets. They found that the mass-ratio 
function of microlensing planets increases as q decreases un- 
til a break at q ~ 1.7 x 1074, below which the planetary 
occurrence rate likely drops. This break suggests that the 
Neptune-mass planets are likely to be the most common of 
cold planets. However, the Suzuki et al. (2016) sample only 
contains two q < 1074 and thus may be affected by small 
number statistics. To examine the existence of the break, a 
larger q < 1074 sample is needed. 

After its commissioning season in 2015, the new-generation 
microlensing survey, the Korea Microlensing Telescope Net- 
work (KMTNet, Kim et al. 2016), has been conducting near- 
continuous, wide-area, high-cadence surveys for ~ 96 deg”. 
The fields with cadences of T > 2 hr + are the KMTNet 
prime fields (~ 12 deg?) and the other fields are the KMT- 
Net sub-prime fields (~ 84 deg”). Since 2016, the detections 
of q < 1074 planets have been greatly increased in two 
ways, and the KMTNet data played a major or decisive role 
in all detections. First, more than ten q < 1074 planets have 
been detected from by-eye searches, including three with 
q < 2 x 1075 (Gould et al. 2020; Yee et al. 2021; Zang 
et al. 2021a). Second, Zang et al. (2021b, 2022a) developed 
the KMTNet AnomalyFinder algorithm to systematically 
search for planetary signals. This algorithm has been applied 
to the 2018 and 2019 KMTNet prime fields (T > 2 hr 3) 
and uncovered five new q < 1074 planets (Zang et al. 2021b; 
Hwang et al. 2022; Gould et al. 2022). Moreover, the system- 
atic search opens a window for a homogeneous large-scale 
KMTNet planetary sample. According to the experience 
from 2018 and 2019 KMTNet prime fields, we expect to de- 
tect > 20 planets with q < 1074 from 2016-2019 seasons. 


! They are OGLE-2005-BLG-169Lb (Gould et al. 2006), OGLE-2005- 
BLG-390Lb (Beaulieu et al. 2006), OGLE-2007-BLG-368Lb (Sumi et al. 
2010), MOA-2009-BLG-266Lb (Muraki et al. 2011), OGLE-2013-BLG- 
0341Lb (Gould et al. 2014b), OGLE-2015-BLG-1670 (Ranc et al. 2019). 


This will be an order of magnitude larger than the Suzuki 
et al. (2016) sample at q < 10-4. 

To build the first KMTNet q < 1074 statistical sam- 
ple, we applied the KMTNet AnomalyFinder algorithm 
to the 2016-2019 KMTNet microlensing events. In this 
paper, we introduce seven new q < 10-^ events from 
this search. They are KMT-2017-BLG-1194, KMT-2017- 
BLG-0428, KMT-2019-BLG-1806/OGLE-2019-BLG-1250, 
KMT-2017-BLG-1003, KMT-2019-BLG-1367, OGLE- 
2017-BLG-1806/KMT-2017-BLG-1021, and KMT-2016- 
BLG-1105. Together with 17 already published and three 
that will be published elsewhere, the KMTNet Anoma- 
lyFinder algorithm found 27 events that can be fit by 
q < 1074 models from 2016-2019 KMTNet data. How- 
ever, whether a planet can be used for statistical studies 
requires further investigations, which is beyond the scope of 
this paper. 

The paper is structured as follows. In Section 2, we briefly 
introduce the KMTNet AnomalyFinder algorithm and the 
procedure to form the q < 1074 sample. In Sections 3, 4 
and 5, we present the observations and the analysis of seven 
q < 10~* events. Finally, we discuss the implications from 
the 2016-2019 KMTNet q < 1074 planetary sample in Sec- 
tion 6. 


2. THE BASIC OF ANOMALYFINDER AND THE 
PROCEDURE 


Section 2 of Zang et al. (2021b) and Section 2 of Zang et al. 
(2022a) together introduced the KMTNet AnomalyFinder 
algorithm. The AnomalyFinder uses a Gould (1996) 2- 
dimensional grid of (to, tef) to search for and fit anoma- 
lies from the residuals to a point-source point-lens (PSPL, 
Paczyński 1986) model. Here to is the time of maximum 
magnification, and teg is the effective timescale. For our 
search, the shortest teg is 0.05 days and the longest teg is 
6.65 days. The parameters that evaluate the significance of 
a candidate anomaly are Ax? and Ax2,,. See Equation (4) 
of Zang et al. (2021b) for their definitions. The criteria of 
Ax? and Ax2,, are the same as the criteria used in Zang 
et al. (2022a); Gould et al. (2022); Jung et al. (2022), with 
Ax? > 200, or Ax? > 120 and Ax2,, > 60 for the KMT- 
Net prime-field events and Ax2 > 100, or Ax2 > 60 and 
Ax2,, > 30 for the KMTNet sub-prime-field events. Future 
statistical studies should use the same criteria. In addition, an 
anomaly is required to contain at least three successive points 
> 2c away from a PSPL model. 

As a result, we found 464 and 608 candidate anomalies 
from 2016-2019 KMTNet prime-field and sub-prime-field 
events, respectively. We checked whether the data from other 
surveys are consistent with the KMTNet-based anomalies 
and cross-checked with C. Han's modeling. We fitted all the 
q < 107° candidates with online data and found 13 new 


candidates with q < 2 x 1074. Then, we conducted tender- 
loving care (TLC) re-reductions and re-fitted the 13 events. 
Of these, eight events unambiguously have q < 1074, three 
events, KMT-2016-BLG-1307, KMT-2017-BLG-0849, and 
KMT-2017-BLG-1057, have 107^ < q < 2 x 107+, and 
two events, KMT-2016-BLG-0625 (Shin et al. in prep) and 
OGLE-2017-BLG-0448/KMT-2017-BLG-0090 (Zhai et al. 
in prep), have ambiguous mass ratios at 1075 < q X 1073 
and will be published elsewhere. 

Among the eight unambiguous q < 1074 events, one 
event, OGLE-2016-BLG-0007/MOA-2016-BLG-088/KMT- 
2016-BLG-1991, will be published elsewhere because it has 
the lowest-q of this sample. We analyze and publish the re- 
maining seven events in this paper. We note that the planetary 
signals of the seven events are not strong, although they are 
confirmed by at least two data sets. We thus further check 
whether the light curves have other similar anomalies, to ex- 
clude the possibility of unknown systematic errors. We ap- 
plied the AnomalyFinder algorithm to the re-reduction data. 
For all of the seven events, besides the known planetary sig- 
nals no anomaly with Ax2 > 20 was detected. Therefore, 
the light curves of the seven events are stable and planetary 
signals are reliable. 


3. OBSERVATIONS AND DATA REDUCTIONS 


Table 1 lists the basic observational information for the 
seven events, including event names, the first discovery 
date, the coordinates in the equatorial and galactic sys- 
tems, and the nominal cadences (T). The seven plane- 
tary events were all identified by the KMTNet post-season 
EventFinder algorithm (Kim et al. 20182). Of them, KMT- 
2019-BLG-1806/OGLE-2019-BLG-1250 and OGLE-2017- 
BLG-1806/KMT-2017-BLG-1021 were discovered by the 
KMTNet alert-finder system (Kim et al. 2018b) and the Early 
Warning System (Udalski et al. 1994; Udalski 2003) of the 
Optical Gravitational Lensing Experiment (OGLE, Udalski 
et al. 2015), respectively, during their observational seasons. 
Hereafter, we designate KMT-2019-BLG-1806/OGLE- 
2019-BLG-1250 and OGLE-2017-BLG-1806/KMT-2017- 
BLG-1021 by their first-discovery name, KMT-2019-BLG- 
1806 and OGLE-2017-BLG-1806. During the 2019 observa- 
tional season, the KMTNet alert-finder system also discov- 
ered KMT-2019-BLG-1367. In addition, OGLE observed the 
locations of KMT-2019-BLG-1367 and KMT-2016-BLG- 
1105 but did not alert them. We also include the OGLE data 
for these two events into the light-curve analysis, for which 
the OGLE data confirm the planetary signals found by the 
KMTNet. MOA did not issue alerts for any of the seven 
events, and there were no follow-up data to the best of our 
knowledge. 

KMTNet conducted observations from three identical 
1.6 m telescopes equipped with 4 deg? cameras in Chile 
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(KMTC), South Africa (KMTS), and Australia (KMTA). 
OGLE took data using an 1.3m telescope with 1.4 deg? field 
of view in Chile. For both surveys, most of the images were 
taken in the J band, and a fraction of V-band images were 
acquired for source color measurements. Each KMTNet V- 
band data point was taken one minute before or after one 
KMTNet I-band data point of the same field. 

The KMTNet and OGLE data used in the light-curve anal- 
ysis were reduced using the custom photometry pipelines 
based on the difference imaging technique (Tomaney & 
Crotts 1996; Alard & Lupton 1998): pySIS (Albrow et al. 
2009, Yang et al. in prep) for the KMTNet data, and Woz- 
niak (2000) for the OGLE data. For each event, the KMTC 
data were additionally reduced using the pyDIA photometry 
pipeline (Albrow 2017) to measure the source color. Ex- 
cept for OGLE-2017-BLG-1806 and KMT-2016-BLG-1105 
whose sources are not located in any OGLE star catalog, the 
I-band magnitudes of the other five events reported in this 
paper have been calibrated to the standard /-band magnitude 
using the OGLE-III star catalog (Szymański et al. 2011). 


4. LIGHT-CURVE ANALYSIS 
4.1. Preamble 


Because all seven events contain short-lived deviations 
from a PSPL model, we first introduce the common methods 
for the light-curve analysis. The PSPL model is described 
by three parameters, £o, uo, and tg, which respectively rep- 
resent the time of lens-source closest approach, the closest 
lens-source projected separation normalized to the angular 
Einstein radius 0g, and the Einstein timescale, 


0 
tg E Om = V KM Are, (1) 


= 1 
Hrel 


where k = a ~ 8.144577, My, is the lens mass, and 
(Trel, Hrel) are the lens-source relative (parallax, proper mo- 
tion). In addition, for each data set ?, we introduce two linear 
parameters, (fs; fp,;), to fit the flux of the source and any 
blend flux, respectively. 

We search for binary-lens single-source (2L1S) models for 
each event. A 2L1S model requires four parameters in addi- 
tion to the PSPL parameters, (s, q, œ, p), which respectively 
denote the planet-host projected separation in units of 0g, the 
planet/host mass ratio, the angle between the source trajec- 
tory and the binary axis, and the angular source radius 6, 
scaled to Og, i.e., p = 0/05. 

Although the final results need detailed numerical analysis, 
some of the 2L1S parameters can be estimated by heuristic 
analysis. A PSPL fit excluding the data points around the 
anomaly can yield the three PSPL parameters, to, uo, and 
tg. If an anomaly occurred at tanom, the corresponding lens- 
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Table 1. Event Names, Alerts, Locations, and Cadences for the six planetary events 


Event Name Alert Date RA 32000 Decl. 32000 £ b T(hr™t) 

KMT-2017-BLG-1194 Post Season  18:17:17.31 —25:19:26.18 +6.63 —4.34 0.4 
KMT-2017-BLG-0428 Post Season  18:05:32.46 —28:29:25.01 +2.59 —3.55 4.0 
KMT-2019-BLG-1806 26 Jul2019 18:02:09.01 —29:24:53.60 +1.41 —3.35 1.0 
OGLE-2019-BLG-1250 0.3 
KMT-2017-BLG-1003 Post Season  17:41:38.76 —24:22:26.18 +3.42 +3.15 1.0 
KMT-2019-BLG-1367 27 Jun 2019 18:09:53.12 —29:45:43.96 +1.93 —4.99 0.4 
OGLE-2017-BLG-1806 14 Oct 2017 17:46:29.58 —24:16:20.17 +4.09 +2.26 0.3 
KMT-2017-BLG-1021 1.0 
KMT-2016-BLG-1105 Post Season  17:45:47.34 —26:15:58.93 +2.30 +1.16 1.0 

source offset, Uanom, and a can be estimated by use a code based on the advanced contour integration code 

(Bozza 2010; Bozza et al. 2018), VBBinaryLensing? to 

" tanom — to E AN compute the 2L1S magnification. For each grid point, we 

Uanom = {/ Uo + tg ; le|- [sin i cor search for the minimum x? by Markov chain Monte Carlo 

(2) (MCMC) x? minimization using the emcee ensemble sam- 


Because the planetary caustics are located at the position of 
[s — 571| ~ Uanom, we obtain 


s4 ~ Yom - —— (3) 


where s = s, and s = s_ correspond to the major-image 
(quadrilateral) and the minor-image (triangular) planetary 
caustics, respectively. For two degenerate solutions with sim- 
ilar q but different s, Ryu et al. (2022) suggested that the ge- 
ometric mean of two solutions satisfies 


Smean — Sr. (4) 


In addition, Zhang et al. (2022) suggested a slightly different 
formalism, and Zhang & Gaudi (2022) provided a theoretical 
treatment of it. For a dip-type planetary signal, Hwang et al. 
(2022) pointed out that the mass ratio can be estimated by 


Ataip V? 
a=( a ) Ê | sin al, (5) 


where Ataip is the duration of the dip, and the accuracy of 
Equation (5) should be at a factor of ~ 2 level. 

To find all the possible 2L1S models, we conduct two- 
phase grid searches for the parameters, (log s, log q, a, p). 
In the first phase, we conduct a sparse grid, which consists 
of 21 values equally spaced between —1.0 < logs < 1.0, 
20 values equally spaced between 0° < a < 360°, 61 val- 
ues equally spaced between —6.0 < logg < 0.0 and five 
values equally spaced between —3.5 < logp < —1.5. We 


pler (Foreman-Mackey et al. 2013), with fixed (log q, log s) 
and free (to, uo, tE, p, a). In the second phase, we conduct a 
denser (log s, log q, œ, p) grid search around each local min- 
imum (e.g., Zang et al. 2022b). Finally, we refine the best-fit 
models by MCMC with all parameters free. 

For degenerate solutions, Yang et al. (2022) suggested that 
the phase-space factors can be used to weight the proba- 
bility of each solution. We follow the procedures of Yang 
et al. (2022) and first calculate the covariance matrix, C, 
of (log s, log q, aœ) from the MCMC chain. Then, the phase- 
space factor is 


p = y det(C). (6) 


Because whether a planet and its individual solutions can be 
used for statistical studies requires further investigations, we 
provide the phase-space factors for the event with multiple 
solutions but do not use them to weight or reject solutions. 

We also investigate whether the inclusion of two high- 
order effects can improve the fit. The first is the microlensing 
parallax effect (Gould 1992, 2000, 2004), which is due to the 
Earth's orbital acceleration around the Sun. We fit it by two 
parameters, 7E, and 7,5, which are the north and east com- 
ponents of the microlensing parallax vector 7rg in equatorial 
coordinates, 


— Trel Herel 


TE Z . 7 
un Hel 0) 


? http;//www.fisica.unisa.it/GravitationAstrophysics/VBBinaryLensing. 
htm 


Table 2. 2L1S Parameters for KMT-2017-BLG-1194 


Parameter A B 
x? /dof 928.0/928 950.6/928 
to (HJD') 7942.66 +0.13 7942.59+0.13 
uo 0.256 0.018 0.246 0.011 
tg (days) 47.0 + 2.5 4T.9 x 1.7 
p(107?) « 2.6 «14 
a (rad) 2.505 + 0.013 2.515 + 0.011 
s 0.8063 + 0.0103 0.8055 + 0.0065 
log q —4.582 + 0.058 | —4.585 + 0.074 


Is OGLE 20.28 + 0.08 20.34 + 0.06 


NOTE—The upper limit on p is 30. 


We also fit the ug > 0 and uo < O0 solutions to consider 
the "ecliptic degeneracy" (Jiang et al. 2004; Poindexter et al. 
2005). For four cases in this paper, the parallax contours take 
the form of elongated ellipses, so we report the constraints on 
the minor axes of the error ellipse, (mg), which is approxi- 
mately parallel with the direction of the Earth's acceleration. 
For the major axes of the parallax contours, 7g, | ^ TEN, 
we only report it when the constraint is useful. 

The second effect is the lens orbital motion effect (Batista 
et al. 2011; Skowron et al. 2011), and we fit it by the pa- 


rameter y = (= d 4) , where ds/dt and da/dt represent 


s ^ dt 
the instantaneous changes in the separation and orientation 
of the two components defined at to, respectively. To exclude 
unbound systems, we restrict the MCMC trials to 6 < 1.0. 
Here f is the absolute value of the ratio of projected kinetic 
to potential energy (An et al. 2002; Dong et al. 2009), 


KE, | _ &Moyr? TE » E à (8) 
PE,| 841? Of "Tg mT78/Ó0g ' 
and where 7g is the source parallax estimated by the mean 


distance to red clump stars in the direction of each event 
(Nataf et al. 2013). 


s=] 


4.2. “Dip” Anomalies 
4.2.1. KMT-2017-BLG-1194 


Figure 1 shows the observed data together with the best-fit 
PSPL and 2L1S models for KMT-2017-BLG-1194. There is 
a dip centered on HJD' ~ 7958.9 (HJD' = HJD—2450000) 
, i.e., tanom ~ 7958.9, with a duration of Atgip ~ 1.05 days. 
The dip and the ridge around the dip are covered by three 
KMTNet sites, so the anomaly is secure. A PSPL fit yields 
(to, uo, tg) = (7942.7, 0.26, 46), and using the heuristic for- 
malism of Section 4.1, we obtain 


a ~ 143.5°; s—5.-0.807; logqg~ —4.68. (9) 


KMT-2017-BLG-1194 


j 1 


Residuals 


191b Vitis -Zs = 135.6 
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Figure 1. The observed data and the 2L1S (the black and orange 
solid lines) and 1L1S models (the grey dashed line) for KMT-2017- 
BLG-1194. The data taken from different data sets are shown with 
different colors. The bottom panels show a close-up of the dip-type 
planetary signal and the residuals to the 2L1S models. 


The grid search yields one solution. Its parameters are 
presented in Table 2 and are in good agreement with the 
heuristic estimates. The top left panel of Figure 2 displays 
the caustic structure and the source trajectory, for which the 
two minor-image planetary caustics are located on both sides 
of the source trajectory. We label the solution as the so- 
lution "A". To further investigate the parameter space and 
check whether the event has the inner/outer solutions (Gaudi 
& Gould 1997), for which the source passes inside (the “In- 
ner" solution) the two planetary caustics (closer to the central 
caustic) or outside (the “Outer” solution), we follow the pro- 
cedures of Hwang et al. (20182). First, we conduct a “hotter” 
MCMC with the error bar inflated by a factor of V3.0. Sec- 
ond, we make a scatter plot of log q versus A£ from the “hot- 
ter" MCMC chain. Here A£ represents the offset between 
the source and the planetary caustic as the source crosses the 
binary axis, 


AE = ugcsc(o) — (s — s 1). (10) 


The resulting scatter plot is shown in Figure 3, from which 
we find another local minimum at A£ ~ 0.02. We label this 
solution as the “B” solution. As shown in the top right panel 
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Figure 2. Geometries of the five “dip” planetary events. In each 
panel, the red lines represent the caustic, the black solid line rep- 
resents the source trajectory, and the line with an arrow indicates 
the direction of the source motion. For the outer solution of KMT- 
2017-BLG-1003, p is constrained at the > 3ø level, so the radius 
of the green dot represents the source radius. For other solutions, p 
only has weak constraints with < 3c, so their source radii are not 
shown. 


of Figure 2, the *B" solution corresponds to the "Inner" so- 
lution. Its parameters from MCMC are given in Table 2 and 
it is disfavored by Ax? = 22.6 compared to the “A” solution. 
In Figure 1, the “B” solution cannot fit the anomaly well and 
all three KMTNet data sets contribute to the Ay?. The ra- 
tio of the phase-space factors is pa : pp = 1 : 0.54, which 
also prefers the "A" solution. Thus, we exclude the “B” so- 
lution. In addition, the models, which have the geometry of 
the “Outer” solution, do not form a local minimum and are 
disfavored by Ax? > 60 compared to the “A” solution. 

For the “A” solution a point-source model is consistent 
within 1c and the 3c upper limit is p « 0.0026. The inclu- 
sion of higher-order effects yields a constraint on Tg, j, and 
with the other 2L1S parameters being almost unchanged. We 
obtain 7g, = —0.18 + 0.35 and adopt the constraints on mg 
and p in the Bayesian analysis of Section 5. This is a new 
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Figure 3. Scatter plot of log q vs. AE for KMT-2017-BLG-1194, 
where A€ = uo csc(o) — (s — s^ 1) represents the offset between 
the source and the center of the planetary caustic at the moment that 
the source crosses the binary axis. The distribution is derived by 
inflating the error bars by a factor of 4/3 and then multiplying the 
resulting x? by 3 for the plot. Red, yellow, magenta, green, blue 
and black colors represent Ax? < 2 x (1,4, 9, 16, 25,00). “A” and 
*B" represent two local minima and the corresponding parameters 
are given in Table 2. 


microlensing planet with q ~ 2.62 x 1075; i.e., about nine 
times the Earth/Sun mass ratio. 


4.2.2. KMT-2017-BLG-0428 


Table 3. 2L1S Parameters for KMT-2017-BLG-0428 


Parameter Inner Outer 
x?/dof 9952.0/9952 9952.1/9952 
to (HJD) 7943.976 + 0.030 7943.978 + 0.031 
uo 0.205 + 0.009 0.205 + 0.009 
tg (days) 44.4+ 1.5 44.3+1.5 
p(10~3) « 6.4 « 6.1 
a (rad) 1.890 + 0.005 1.889 + 0.005 
E 0.8819 + 0.0044 0.9146 + 0.0050 
log q —4.295 + 0.072 —4.302 + 0.075 
Is,0GLE 20.43 + 0.05 20.43 + 0.05 


Figure 4 shows a AI ~ 0.12 mag dip at tanom ~ 7947.00, 
with a duration of Atgi, ~ 0.74 days. The dip is defined by 
the KMTA and KMTC data, and the subtle ridges are sup- 
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Figure 4. The observed data and models for KMT-2017-BLG-0428. 
The symbols are similar to those in Figure 1. In the top panel, the 
black arrow indicates the position of the planetary signal. 


ported by both the KMTC and KMTS data. These data were 
taken in good seeing (1."4—2."5) and the anomaly does not 
correlate with seeing, sky background or airmass. In addi- 
tion, Ishitani Silva et al. (2022) found that the KMTA data 
show systematic errors and excluded them from the analysis. 
In that case, the KMTA data exhibit similar residuals from 
one-night data in many places of the light curves. For the 
present case, the anomaly is mainly covered by the KMTA 
data, but as presented in Section 2, there is no similar devi- 
ation in other places of the light curves. We also carefully 
checked the KMTA data but did not find any similar resid- 
uals. Hence, the anomaly is secure. Applying the heuristic 
formalism of Section 4.1, we obtain 

a ~ 108.4^; s=s_~0.898; logq~ —4.19. (11) 

The 2L1S modeling yields two degenerate solutions with 
Ax? = 0.1. As shown in Figure 2, the two solutions are sub- 
jected to the inner/outer degeneracy. Their parameters are 
given in Table 3, for which a and q are consistent with Equa- 
tion (11). For s, we note that the geometric mean of the two 
solutions, Smean = 0.898 + 0.005, is in good agreement with 
Equation (11) and thus the formalism of Ryu et al. (2022). 
In addition, the observed data only provide a 3c upper limit 
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on p, and a point-source model is consistent within 1c. The 
ratio of the phase-space factors iS Pinner : Pouter = 0.78 : 1. 

With high-order effects, we find that the x? improvement is 
~ 3 and other parameters are almost the same. The constraint 
of TE, Tg | = —0.35 + 0.26, will be used in the Bayesian 
analysis. This is a microlensing planet with a Neptune/Sun 
mass ratio. 


4.2.3. KMT-2019-BLG-1806 
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Figure 5. The observed data and models for KMT-2019-BLG- 1806. 
The symbols are similar to those in Figure 1. In the top panel, the 
black arrow indicates the position of the planetary signal. 


The anomaly of KMT-2019-BLG-1806 is also a dip, as 
shown in Figure 5. The dip has Atg;, ~ 0.6 days and cen- 
ters on tanom ~ 8717.72. The dip is defined by the KMTC 
data and the two contemporaneous OGLE points, which were 
taken in good seeing (1."1—2."4) and low sky background. 
Hence, the anomaly is secure. Applying the heuristic for- 
malism of Section 4.1, we obtain 


a ~ 123.3°; s—5.7-0.985; logq ~ —4.56. (12) 


In addition, given the Einstein timescale (tg ~ 135 days), 
we expect that 7g should be either measured or strongly con- 
strained. 

The 2L1S modeling also finds a pair of inner/outer solu- 
tions and combined the ug > 0 and uo < 0 degeneracy 
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Table 4. 2L1S Parameters KMT-2019-BLG- 1806 


Parameter Inner Outer 
uo > 0 uo <0 uo > O uo <0 
x?/dof 3132.5/3132 3132.9/3132 3132.2/3132 3131.8/3132 
to (HJD) 8715.452 + 0.015 8715.451 +0.015 8715.453+0.014 8715.453 + 0.015 
uo 0.0260 + 0.0017 —0.0251 + 0.0020 0.0257 + 0.0016 —0.0255 + 0.0015 
tg (days) 132.8 + 8.1 138.5 + 10.8 134.1 + 7.9 135.6 + 7.9 
p(1073) «18 «1.8 «1.9 <1.7 
a (rad) 2.150 + 0.008 —2.147 + 0.008 2.151 + 0.009 —2.148 + 0.008 
E 0.9377 + 0.0069 0.9383 + 0.0073 1.0339 + 0.0084 1.0352 + 0.0085 
log q —4.724 + 0.117 —4.734 + 0.109 —4.717 + 0.117 —4.714 + 0.116 
TE,N —0.055 + 0.150 —0.066 + 0.161 —0.060 + 0.156 —0.019 + 0.160 
TE,E —0.058 + 0.017 —0.059 + 0.014 —0.057 + 0.017 —0.060 + 0.013 
Is 21.33 + 0.07 21.37 + 0.09 21.34 + 0.07 21.35 + 0.07 


T 
KMTA19 ] 


there are four solutions in total. See Table 4 for their param- 
eters. The inclusion of mg improves the fits by Ax? = 31, 
and all four data sets contribute to the improvement, so the 172 
parallax signal is reliable. The angle of the minor axis of 
the parallax ellipse (north through east) is i» = 82.0? and 
i = 82.5? for the uo > 0 and ug < O solutions, respec- 
tively. mp = 0.06 + 0.01, and zg,, is constrained to UMP 
be o(7p,1) ~ 0.2. We obtain Smean = 0.985 + 0.008, isopE HT T 
a = 123.1 + 0.5, and logq ^ —4.72, in good agreement i 
with Equation (12). The ratio of the phase-space factors is E 
Dinner : Pouter = 0.80 : 1. g 

We find that the inclusion of the lens orbital motion effect 
only improves the fit by Ax? < 1 for 2 degree-of-freedom 
and is not correlated with mg, so we exclude the lens orbital 
motion effect. With q ~ 1.9 x 1075, this new planet is the 
fifth robust q < 2 x 10^? microlensing planet. 
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4.2.4. KMT-2017-BLG-1003 


Mis Xs = 247.8 


Figure 6 shows the light curve and the best-fit models for 
KMT-2017-BLG-1003. The KMTC data show a sudden dip 
and the ridge is confirmed by the KMTC and KMTS data. 
These data were taken in good seeing (1."2-2."2) and low 
sky background, so the anomaly is of astrophysical origin. 
Although the end of the dip is not covered, the KMTC point 
at HJD’ = 7870.66 indicates Ataip < 0.85 days, which 
yields 
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Figure 6. Light curve and models for KMT-2017-BLG-1003. The 
symbols are similar to those in Figure 1. 


0761.22; s=s_~0.903; logq < —3.6. (13) 


The numerical analysis yields a pair of inner/outer solu- 
tions, and Table 5 lists their parameters. As shown in Fig- 
ure 2, the “Outer” solution has caustic crossings, so its p is 
measured at the 4.5c level. For the “Inner” solution, a point- 
source model is consistent within 2c. We note that the geo- 
metric mean of s, Smean = 0.899 + 0.004, which is slightly 


different from s_ by 1c. This indicates that the prediction of 
Ryu et al. (2022) might be imperfect for minor-image anoma- 
lies with finite-source effects or incomplete coverage. The 
ratio of the phase-space factors is Pinner : Pouter = 0.80 : 1. 


Table 5. 2L1S Parameters for KMT-2017-BLG-1003 


Parameter Inner Outer 
x7 /dof 2433.2/2433 2433.0/2433 
to (HJD’) 7872.484 + 0.020 7872.482 + 0.020 
uo 0.179 0.005 0.179 + 0.005 
te (days) 25.65 + 0.57 25.66 + 0.59 
p(107?) « 6.7 5.22 + 1.16 
a (rad) 1.073 + 0.006 1.072 + 0.006 
s 0.8889 + 0.0043 0.9096 + 0.0045 
logq —4.260 + 0.152 — —4.373 + 0.144 
Is,oGLE 19.30 + 0.04 19.30 + 0.04 


With high-order effects, the x? improvement is 1.7. Al- 
though this event is shorter than the first two events, mg is 
better constrained due to the about one magnitude brighter 
data, with TE, = —0.11-£0.15. This is another Neptune/Sun 
mass-ratio planet. 


4.2.5. KMT-2019-BLG-1367 
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Figure 7. Light curve and models for KMT-2019-BLG-1367. The 
symbols are similar to those in Figure 1. 


Table 6. 2L1S Parameters for KMT-2019-BLG-1367 


Parameter Inner Outer 
x7 /dof 1404.0/1404 1404.2/1404 
to (HJD) 8667.883 + 0.051 8667.884 + 0.048 
uo 0.083 0.009 0.082 0.009 
tp (days) 39.3 + 3.8 39.8 + 4.0 
p107?) < 5.3 < 5.6 
a (rad) 1.208 + 0.016 1.207 + 0.016 
E 0.9389 + 0.0066 0.9763 + 0.0070 
log q —4.303 + 0.118 —4.298 + 0.103 
Is,0GLE 21.46 + 0.13 21.48 + 0.13 


Figure 7 shows a dip 1.2 days before the peak of an oth- 
erwise normal PSPL event, with a duration of Atgip ~ 0.35 
days. The dip-type anomaly is covered by the KMTC data 
and one contemporaneous OGLE point, and these data were 
taken in good seeing (< 2."0) and low sky background. 
Therefore, the anomaly is secure. Applying the heuristic for- 
malism of Section 4.1, we obtain 


a ~ 69.5°; s=s_~0.957; logq ~ —4.33. (14) 


The 2L1S modeling also yields a pair of inner/outer solu- 
tions, with Ay? = 0.2. The resulting solutions are given 
in Table 6 and Figure 2. A point-source model is consistent 
within 1c and the 3c upper limit is p < 0.0056, so we expect 
that the Ryu et al. (2022) formula is applicable. We obtain 
Smean = 0.957 + 0.007, in good agreement with s... The 
ratio of the phase-space factors is Dinner : Pouter = 0.82 : 1. 
We find that the inclusion of higher-order effects only im- 
proves the fitting by Ay? < 1 and the 1c uncertainty of 
parallax is > 0.9 at all directions, so the constraint on mg is 
not useful for the Bayesian analysis. This is another planet 
with a Neptune/Sun mass ratio. 


4.3. "Bump" Anomalies 


For bump-type planetary signals, we also check whether 
the observed data can be fitted by a single-lens binary-source 
(1L2S) model (Gaudi 1998) because it can also produce such 
anomalies (e.g., Hwang et al. 2013; Jung et al. 2017; Rota 
et al. 2021). For a 1L2S model, its magnification, A», is the 
superposition of magnifications for two single-lens single- 
source (1L1S) models, 

A Ai fia + Aofax — Ai +qfaA2. _ far 
A = , qf = , 
fux + Jaa l-qgf JLA 

(45) 


where f; is the source flux at wavelength A, and 7 = 1 and 
i — 2 correspond to the primary and the secondary sources, 
respectively. 
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Figure 8. Light curve and models for OGLE-2017-BLG-1806. The 
symbols are similar to those in Figure 1. Different with the previous 
four events, the anomaly is bump-type, so the best-fit 1L2S model 
is provided. 


As shown in Figure 8, the light curve of OGLE-2017-BLG- 
1806 exhibits a bump centered on tanom ~ 8003.5, defined 
by the KMTC and KMTS data. Except for two KMTS points, 
all the KMTC and KMTS data during 8003 < HJD' « 8005 
were taken in good seeing (« 2."2) and low sky back- 
ground. In addition, most of the data before the bump 
(8000 < HJD' « 8003) are fainter than the 1L1S model. 
Hence, the signal is secure. Because both the major-image 
and the two minor-image planetary caustics can produce a 
bump-type anomaly (e.g., Wang et al. 2022), we obtain 

a ~ 184.6? for s, ~ 1.15; a~ 4.6? for s_ ~ 0.86. 

(16) 

The grid search returns three local minima, and their caus- 
tic structures are given in Figure 9. As expected, the three so- 
lutions respectively correspond to sources crossing a major- 
image (quadrilateral) planetary caustic and two minor-image 
(triangular) planetary caustics. We label the three solutions 
as “Close A", “Close B", and “Wide”, respectively, and their 
parameters are presented in Table 7. 
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Figure 9. Geometries of OGLE-2017-BLG-1806. The symbols are 
similar to those in Figure 2. For the two “Close” solutions, p is 
constrained at the > 3ø level, so the radius of the two green dots 
represent the source radius. For the “Wide” solution, p only has 
weak constraints with < 30, so its source radius is not shown. 


The “Close A" solution provides the best fit to the ob- 
served data, and the “Close B" and “Wide” solutions are 
disfavored by Ax? — 14.1 and 8.3, respectively. We find 
that the inclusion of the parallax effect improves the fit by 
Ax? = 7.8, 20.4, and 11.1 for the “Close A", “Close B", 
and “Wide” solutions, respectively, and during the anomaly 
region (7998 < HJD' « 8008), Ax? =2.2, 22.3, and 6.8. 
With the anomaly removed, fitting the data by a 1L1S model 
yields a similar constraint on Tg, and a weaker constraint on 
TE, 1, With a (75,1) ~ 0.5. Thus, the long planetary signal 
plays an important role in improving the constraint on 7g, 1 
and reduces the x? differences between the three solutions. 

The ratio of the phase-space factors is pcioseA : PCloseB : 
Dwide = 1 : 0.95 : 0.61. For the “Close A", and “Close B" 
solutions, the bump was produced by a caustic crossing, so o 
is constrained at the 3c level. For the “Wide” solution, the 
bump was a result of a cusp approach. Although the “Wide” 
solution has caustic crossing features, due to the lack of data 
during the crossing, a point-source model is consistent within 
lo. 

The 1L2S model is disfavored by Ay? = 30.7 compared 
to the “Close A” solution, and the 1L2S parameters are shown 
in Table 8. Although the 1L2S model fits the bump well, it 
provides a worse fit to the observed data before the bump, 
during which most of the data from the three KMTNet sites 


Table 7. 2L1S Parameters for OGLE-2017-BLG-1806 
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Parameter Close A Close B Wide 
uo > 0 uo <0 uo > 0 uo < 0 uo > 0 uo «O0 
x?/dof 1650.9/1651 1650.7/1651 1664.8/1651 1665.5/1651 1659.1/1651 1659.0/1651 
to (HJD') 8024.392 + 0.020 8024.393 + 0.019 8024.388 + 0.020 8024.388 + 0.020 8024.388 + 0.020 8024.379 + 0.020 
uo 0.0249 + 0.0016 — —0.0260 + 0.0016 0.0256 + 0.0020  —0.0253 + 0.0019 0.0269 + 0.0018 —0.0257 + 0.0017 
tg (days) 69.4 + 4.0 66.8 + 3.9 69.4 + 4.8 69.6 + 4.6 64.5 + 3.9 67.0 + 3.9 
p(1073) 1.745075 1.83* 0:20 1.50736? 1.657087 « 2.8 «24 
a (rad) 0.001 + 0.034 —0.002 + 0.037 0.267 + 0.066 —0.263 + 0.068 3.121 + 0.034 —3.121 + 0.036 
E 0.8609 + 0.0069 0.8566 + 0.0075 0.8592 + 0.0085 0.8601 + 0.0080 1.1900 + 0.0117 1.1806 + 0.0108 
log q —4.392 + 0.180 —4.352 + 0.171 —4.766 + 0.220 —A4.768 + 0.209 —4.317 + 0.126 —4.44] + 0.168 
TE,N —0.278 + 0.148 0.292 + 0.170 0.774 + 0.315 —0.756 + 0.326 —0.535 + 0.175 0.504 + 0.170 
TE,E 0.105 + 0.056 0.144 + 0.058 0.157 + 0.070 0.124 + 0.059 0.120 + 0.065 0.133 + 0.056 
Is KMTC 21.12 + 0.07 21.07 + 0.07 21.10 + 0.08 21.10 + 0.08 21.03 + 0.07 21.08 + 0.07 
Table 8. 1L2S Parameters for OGLE-2017-BLG-1806 and KMT-2016-BLG-1105 
Parameters OGLE-2017-BLG-1806 KMT-2016-BLG-1105 
uo > 0 uo < 0 
x?/dof 1682.0/1651 1681.4/1651 2298.7 /2288 
to,1 (HJD’) | 8024.383 + 0.020 8024.381 + 0.020 7555.972 + 0.094 
to,2 (HJD’) | 8003.876 + 0.274 8003.913 + 0.253 7547.890 + 0.021 
Uo,1 0.0288 + 0.0023 — —0.0282 + 0.0019 0.143 + 0.022 
0,2 0.003 + 0.025 —0.004 + 0.023 0.0001 + 0.0007 
tg (days) 61.2 + 4.3 62.2 + 3.5 44.9 d 5.8 
p2 (107?) Eds, < 7.0 < 3.3 
qr. 1(107?) 2.76 + 0.76 2.63 + 0.74 1.98 + 0.48 
TEN 0.041 + 0.388 0.059 + 0.355 
TEE 0.111 + 0.072 0.117 + 0.063 
Is,kMTC 20.96 + 0.09 20.98 + 0.07 21.31 + 0.18 


are fainter than the 1L2S model. Hence, the 1L2S model is 
rejected. We find that the lens orbital motion effect is not 
detectable (Ax? < 0.5), so we adopt the parameters with the 
microlensing parallax effect as our final results. 


4.3.2. KMT-2016-BLG-1105 


The anomaly in Figure 10 is a short-lived bump centered 
On tanom ^ 7547.85, which is defined by four KMTC data 
points and supported by one OGLE data point. These data 
were taken in good seeing (« 2."0) and low sky background, 
so the anomaly is secure. Similar to OGLE-2017-BLG-1806, 
we expect both the major-image and the minor-image plane- 


tary caustics can produce the bump and obtain 


a ~ 39.3? for s_ ~ 0.89. 
(17) 
The 2L1S modeling yields five solutions, including one 
with the minor-image planetary caustics and four with the 
major-image planetary caustics. We label them as "Close", 
“Wide A’, “Wide B", “Wide C" and “Wide D", respectively, 
and their parameters are given in Table 9. Figure 11 displays 
the caustic structures and source trajectories. The “Wide 
A’, “Wide B” and “Close” solutions exhibit caustic cross- 
ings, but only for the “Wide B” and “Close” solutions p are 
constrained at the > 3o level. For the “Wide A", “Wide C" 
and “Wide D” solutions, a point-source model is consistent 


a ~ 219.3? for s} ~ 1.13; 
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Table 9. 2L1S Parameters for KMT-2016-BLG-1105 


Parameters Wide A Wide B Wide C Wide D Close 
x?Idof 2286.7/2288 2289.0/2288 2291.1/2288 2289.4/2288 2290.2/2288 
to (HJD/) | 7555.834+ 0.096 7555.789 0.102 7555.772 -- 0.093 7555.781 20.099 7555.896 + 0.093 

uo 0.171 + 0.012 0.153 + 0.013 0.154 + 0.014 0.154 + 0.014 0.148 + 0.008 
tg (days) 38.8 + 2.0 42.4 + 2.9 42.5 3: 3.1 42.4 3: 3.1 43.3 1.8 
pı (107°) < 2.4 2.92 + 0.82 < 4.6 < 5.5 0.75 + 0.14 
a (rad) 3.836 + 0.014 3.830 + 0.016 3.832 + 0.014 3.831 + 0.014 0.691 + 0.021 
E 1.143 + 0.009 1.136 + 0.011 1.155 + 0.012 1.106 + 0.013 0.888 + 0.007 
log q —5.194 + 0.248 —4.423 + 0.197 —4.069 + 0.182 —4.184 + 0.206 —5.027 + 0.080 
Is kmTC 21.09 + 0.08 21.20 + 0.05 21.22 + 0.11 21.22 £0.11 21.27 + 0.06 
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Figure 10. Light curve and models for KMT-2016-BLG-1105. The 
symbols are similar to those in Figure 1. Because a 1L2S model can 
produce a short-lived bump, the best-fit 1L2S model is also shown. 


within Ay? = 3, 1, and 1, respectively, and thus we only 
report their 3c upper limit on p in Table 9. The ratio of the 
phase-space factors is pwideA : PWideB : DWideC : PWideD : 
DClose — 0.82 : 0.76 : 0.74 : 1 : 0.41, so the wide solutions 
are slightly favored in the phase-space factors. 

For the "Close" solution, the bump was produced by a 
cusp approach with the lower triangular planetary caustic, 
followed by a dip that occurred in the data gap between 
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Figure 11. Geometries of KMT-2016-BLG-1105. The symbols are 
similar to those in Figure 2. 


HJD' — 7548.0 and HJD' — 7548.3. If the bump were 


produced by a cusp approach with the upper triangular plan- 
etary caustic, there would be a dip before the bump, but the 
region before the bump is well covered by the KMTS and 
the KMTC data, which are consistent with the 1L1S model. 
Thus, the minor-image perturbation only has one solution. 
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Figure 12. Scatter plot of log q vs. A£ for KMT-2016-BLG-1105. 
The distribution is derived by inflating the error bars by a factor of 
v2.5 and then multiplying the resulting x? by 2.5 for the plot. The 
colors are the same as those in Figure 3. "A", “B”, "C", and “D” 
represent four local minima and the corresponding parameters are 
given in Table 9. 


For the four “Wide” solutions, the “Wide A" and “Wide B" 
solutions have a source crossing the planetary caustic, and 
the “Wide C" and “Wide D" solutions that contain a source 
that passes to one side or the other of the planetary caus- 
tic. This topology is qualitatively similar to the topology of 
OGLE-2017-BLG-0173 (Hwang et al. 20182). We thus also 
investigate the parameter space by a "hotter" MCMC with the 
error bar inflated by a factor of \/2.5. The resulting scatter 
plot is shown in Figure 12, from which we find that the topol- 
ogy of KMT-2016-BLG-1105 has differences in three aspects 
from the topology of OGLE-2017-BLG-0173. First, for the 
two solutions in which the source passes to one side or the 
other of the planetary caustic, OGLE-2017-BLG-0173 has 
caustic crossings and the source is comparable to the size of 
the planetary caustic, but in the present case, the source does 
not cross the caustic. Second, for the solution in which the 
source passes directly over the planetary caustic, the source is 
much larger than the planetary caustics in the case of OGLE- 
2017-BLG-0173, while the source of KMT-2016-BLG-1105 
is smaller than the caustic. Third, OGLE-2017-BLG-0173 
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exhibits a bimodal minimum when the source passes directly 
over the caustic, and the mass-ratio difference between the 
two local minima is Alogq « 0.1. The corresponding solu- 
tions for KMT-2016-BLG-1105, the “Wide A" and “Wide B" 
solutions, have Alogq ~ 1. We note that the “Wide A" and 
“Wide B” solutions have AE ~ 0.00 and —0.01, respectively. 
Considering the approximate symmetry with respect to A£, 
one might expect an additional minimum that has A£ ~ 0.01 
and a similar logg as the logq of the “Wide B" solution. 
However, such a potential solution “disappeared” from the 
numerical analysis. Because the trajectories of the “Wide” B 
solution and the putative minimum at A£ ~ 0.01 should be 
almost symmetric with respect to the center of the caustics, 
their corresponding planetary signals should also be almost 
symmetric. As shown in Figure 10, the “Wide B" solution 
drops rapidly during the caustic exit, followed by a dip, so 
the putative minimum at A£ ^ 0.01 should contain a dip fol- 
lowed by a sudden rise during the caustic entry, which is not 
supported by the KMTC and KMTS data. Thus, in Figure 12 
this topology is absorbed into the MCMC chain of the “Wide 
D" solution and there is no new discrete solution. 

We also check whether the bump-type anomaly can be 
fitted by a 1L2S model. Table 8 lists the IL2S parame- 
ters. We find that the best-fit 1L2S model is disfavored by 
Ax? = 12.0 compared to the best-fit 2L1S model. The best- 
fit IL2S model has p; = 0.0018. We note that the flux ra- 
tio is qf, ~ 2 x 107%, corresponding to a magnitude dif- 
ference of 6.7 mag. According to Section 5, the primary 
source lies 4.1 mag below the red giant clump, so the pu- 
tative source companion would have an absolute magnitude 
of Mr.» ~ 10.7 mag, corresponding to an angular source ra- 
dius of 042 ~ 0.1 pas. This yields the lens-source relative 
proper motion of fire = 0. 2/p2/tg ~ 0.5 mas yr. !, which 
is lower than the typical Hre] of bulge microlensing events 
(See Figure 2 of Zhu et al. (2017) for examples). However, 
a model with o» = 0 is only disfavored by Ax? = 1, so any 
reasonable 144; is only disfavored by Ax? < 1. Thus, while 
the planetary model is strongly favored, there is a possibility 
that the anomaly is caused by a second source. 

With high-order effects, we find that Ax? < 1 and the 
lo uncertainty of parallax is > 0.9 at all directions, so the 
constraint on 7rg is not useful. 


5. SOURCE AND LENS PROPERTIES 
5.1. Preamble 


Combining Equations (1) and (7), the mass My, and dis- 
tance Dy of a lens system are related to the angular Einstein 
radius 0g and the microlensing parallax 7g by (Gould 1992, 
2000) 

On . au 


Dy, = ————— 1 
KTE H Tele + Tg’ (18) 
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Figure 13. Color magnitude diagrams for the seven planetary events analyzed in this paper. The first five CMDs are constructed using the 
OGLE-III star catalog (Szymański et al. 2011), and the other two CMDs are constructed using the KMTC pyDIA photometry reduction. For 
each panel, the red asterisk and the blue dot are shown as the centroid of the red giant clump and the microlensed source star, respectively. The 
three green dots on the CMDs of KMT-2017-BLG-1194, KMT-2019-BLG-1806, and KMT-2017-BLG-1003 represent the blended light. For 
the bottom panel, the yellow dots represent the HST CMD of Holtzman et al. (1998) whose red-clump centroid has been matched to that of 
KMTC using (V — J, I)a,usr = (1.62, 15.15) (Bennett et al. 2008). 


Table 10. CMD parameters, 0.., 0g and rei for the five “dip” planetary events 


Parameter KB171194 KB170428 KB191806 KB171003 KB191367 
Inner Outer 
(V -I)a 1.82 + 0.01 1.95 + 0.01 2.23 + 0.01 2.39 + 0.01 + 1.70 + 0.01 
la 15.25 +0.01 | 15.39-20.01 | 15.7923: 0.02 | 16.04 0.01 1} 15.13 + 0.01 
Ia,o 14.26 + 0.04 | 14.36 + 0.04 | 14.39 +0.04 | 14.34 c 0.04 + 14.37 + 0.04 
(V — I)s 1.47 + 0.07 1.95 + 0.04 1.93 + 0.03 2.00 + 0.02 + 1.70 + 0.03 
Is 20.28 + 0.08 | 20.43 +0.05 | 21.35 +0.07 | 19.30 +0.04 19.30 +0.04 | 21.47 + 0.13 
(V — I)s,o 0.71 + 0.08 1.06 + 0.05 0.76 + 0.05 0.67 + 0.04 + 1.06 + 0.04 
Is,0 19.29+0.09 | 19.40+0.07 | 19.95+0.08 | 17.60+0.06 17.60+0.06 | 20.71 0.14 
0. (uas) 0.448 + 0.038 | 0.578 + 0.034 | 0.345 + 0.020 | 0.942 + 0.046 0.942 + 0.046 | 0.316 + 0.023 
0g (mas) > 0.17 > 0.09 > 0.19 > 0.14 0.180 + 0.041 > 0.06 
[re (mas yr!) > 1.3 > 0.74 > 0.51 > 2.0 2.56 + 0.58 > 0.53 


NOTE—(V — I)a,o = 1.064 


Table 11. CMD parameters, 0., 0g and Hre] for OGLE-2017-BLG-1806. 


t 0.03. Event names are abbreviations, e.g., KMT-2017-BLG-1194 to KB171194. 


Parameter Close A Close B Wide 
uo > 0 uo < 0 uo > 0 uo < 0 uo > 0 uo < 0 
(V — Da 2.89 + 0.01 c— c— + + c— 
la 16.42 + 0.02 + + + + + 
Ico 14.33 + 0.04 + +} + + 4 
(V — I)s 2.66 + 0.03 c— c c— c— c— 
"n 21.12+0.07 21.07 +0.07 | 21.10 +0.08 21.10 +0.08 | 21.03 +0.07 21.08 + 0.07 
(V — D)s,o 0.83 + 0.04 + + + — — 
Iso 19.03 + 0.08 18.98 + 0.08 19.01 + 0.09 19.01 + 0.09 18.94 +0.08 18.99 + 0.08 
6. (uas) 0.561 + 0.031 0.574 0.031 | 0.566 + 0.033 0.566 + 0.033 | 0.584 + 0.032 0.571 + 0.032 
0r (mas) (30702 ah | O37702 — 00349] 00 > 0.21 > 0.24 
Ux (masyr !)) |  1.69*022 Li 1.98+0:65 2.08 10:95 >1.2 >13 
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Table 12. CMD parameters, 0.., 0g and jure: for KMT-2016-BLG-1105. 


Wide A Wide B Wide C Wide D Close 
la 17.20 + 0.01 c — com — 
lao 14.39 + 0.04 1 c 1 + 
Is 21.09 40.08 21.20+0.05 21.22-40.11 21.22-40.11 21.27 X 0.06 
(V — I)s,o 0.74 + 0.07 0.75 + 0.07 0.75 + 0.07 0.75 + 0.07 0.75 + 0.07 
Is, 18.28 +0.09 18.39+0.07 18.41+0.12 18.41 +0.12 18.46 + 0.07 
Ox (uas) 0.732 + 0.057 0.702 + 0.051 0.696 + 0.061 0.696 + 0.061 0.680 + 0.050 
0r (mas) > 0.31 0.240 + 0.070 > 0.15 > 0.13 0.907 + 0.182 
[rel (mas yr“) > 2.9 2.07 + 0.62 > 1.3 >11 7.65 + 1.54 


To obtain the angular Einstein radius through 0g = 0../p, 
we first estimate the angular source radius 0, by locating 
the source on a color-magnitude diagram (CMD, Yoo et al. 
2004). For each event, we construct a V — I versus 1 CMD 
using the ambient stars of the OGLE-III catalog (Szymański 
et al. 2011) or the KMTC images with the pyDIA reductions. 
See Figure 13 for the CMDs of the seven planetary events. 
We estimate the centroid of the red giant clump as (V —I, I)a 
from CMDs and adopt the de-reddened color and magnitude 
of the red giant clump, (V — I, I)a,o, from Bensby et al. 
(2013) and Table 1 of Nataf et al. (2013). We obtain the 
source apparent magnitude from the light-curve analysis of 
Section 4, and the source color by a regression of the KMTC 
V versus I flux with the change of the lensing magnification. 

We find that the V-band observations of KMT-2016-BLG- 
1105 have insufficient signal-to-noise ratio to determine the 
source color, so we estimate the source color by the Hub- 
ble Space Telescope (HST) CMD of Holtzman et al. (1998) 
(see Section 5.7 for details). Finally, using the color/surface- 
brightness relation of Adams et al. (2018), we obtain the an- 
gular source radius 0... Tables 10, 11 and 12 present the CMD 
values and (0., 0g, 1,41) from the procedures above. 

Because the blended light could provide additional con- 
straints on the lens properties (e.g., the upper limits for the 
lens brightness), we also check the brightness and the as- 
trometric alignment of the baseline object. For KMT-2017- 
BLG-0428 and KMT-2019-BLG-1806, we adopt the i’-band 
baseline images taken by the 3.6m Canada-France-Hawaii 
Telescope (CFHT) from 2020 to 2022, whose seeing FWHM 
is 0."55—-0."70. For the other five events which do not have 
any CFHT image, we check the baseline objects from the 
KMTC pyDIA reduction, whose seeing FWHM is about 
1.0. 

Because none of the seven planetary events have simulta- 
neous measurements of mg and Og at the > 3o level, the lens 
masses and distances cannot be directly determined by Equa- 


tion (18). We conduct a Bayesian analysis using a Galactic 
model to estimate the lens properties. The Galactic model 
and the procedures we adopt are the same as described in 
Zang et al. (2021b). We refer the reader to that work for de- 
tails. The only exception is that we include upper limits of 
the lens light, TL limit, from the analysis of the blended light. 
We adopt the mass-luminosity relation of Wang et al. (2018), 
My 
Mr = 4.4 — 8.5 log Mo’ (19) 
where Mpy is the absolute magnitude in the J band, and we 
reject trial events for which the lens properties obey 
Di 
Mr+5 log — + ALD, < Ty, limit; (20) 
10pc 

where Ar. p, is the extinction at D. We adopt an extinction 
curve with a scale height of 120 pc. For the five events with 
OGLE CMDs, the total extinction is derived from the CMD 
analysis, Ay = Ia — Icio. For the other two events with 
KMTC CMDs, we adopt the extinction in the K band from 
Gonzalez et al. (2012) and Ay = 7.26 Ax from Nataf et al. 

(2016). 

Table 13 presents the resulting Bayesian estimates of the 
host mass Mj, ost, the planet mass Mpjanet, the lens distance 
Dy, the projected planet-host separation a, and the lens- 
source relative proper motion Jo. For events with multi- 
ple solutions, we show the results for each solution and the 
"combined results" of combining all solutions weighted by 
their Galactic-model likelihood and exp(—Ay?/2), where 
Ax? is the y? difference compared to the best-fit solu- 
tion. Here the Galactic-model likelihood represents the total 
weight for the simulated events given the error distributions 
of tg, 0g and mg. See Equation (16) of Zang et al. (2021b) 
for the weight procedures. 

We do not adopt the “combined results" as the final phys- 
ical parameters but just show them for consideration, be- 


Table 13. Physical parameters of the six planetary events from a Bayesian analysis. 


Event Solution Physical Properties Relative Weights 
Myost[Mo] Mpianet|Me] Dxikpe] alau]  preif[masyr~'] | Gal.Mod. x? 
KB171194 divo Soa. apo OLET 4.294250 
KB170428 Inner Ost 2- Xo 54035 178105 3271138 0.99 1.00 
Outer Gao. Sots, 5.40133 1:8519% — S983 1.00 0.95 
Combined Cac  5.59t357 ba? 1.814955 3.27122 
KB191806 Inner (uo > 0) | 0.757334 cert. 6621972 2.87109% Lig 1.00 0.70 
Inner (uo < 0) | 0.747332 4.471139 — pg NATUS jd 0.84 0.58 
Outer (uo > 0) | 0.737022 Aon, 6.68972 Gto s Lie 0.98 0.82 
Outer (uo < 0) | 0.754954 ATIE ls, SITO Ags e 0.98 1.00 
Combined Oar? 4.671136 don. 309008 lagu. 
KB171003 Inner te ME vou ET AR o EE us 1.00 0.90 
Outer ov sy mies d 7.5 — EOD 0.74 1.00 
Combined 0.327031 5.191333 7.091981 dq38 3.114035 
KB191367 Inner DT. 4.061358 — AuoR Sao 1.674949 3.92 *257 1.00 1.00 
Outer 0.351915 4.121259 aT. 1731951 eo 0.96 0.90 
Combined 0.254918 — pato 4.671245 1.704959 sor 
OB171806 Close A (uo > 0) | 0.44* 033 Darren Bol. L8 n dto 0.85 0.90 
Close A (uo < 0) | 0.334333 4.831333 6.174922 4597 0 2.131036 1.00 1.00 
Close B (uo > 0) | 0.4470°32 zd as deb tons XU IT 2.997075 0.21 ipe 
Close B (uo < 0) | 0.48703? 2.68119 — ste olt s 2.097095 016 10-3? 
Wide (uo > 0) | 0.34*071 nag (300715. 315975 4.924175 i Joms 
Wide (uo < 0) | 0.417023 Aer LETI, 2 5.48 1-20 024 10°18 
Combined 0.987525 ba dp Le 2.05+9:32 
KB161105 Wide A 0.431922 — 95804 “Set 2.931989 6.481298 0.37 1.00 
Wide B | AGE qusc B08, ^ a ee 1.00 0.32 
Wide C OAs) — qud" 49075 Apr 3.971288 0.66 0.11 
Wide D (.44* 0-27 9.51130 quat Out 4.111285 0.64 0.26 
Close 0.431913 ie E e Xu EE o^ 0.29 0.17 
Combined 0.414925 250", - BOB oat 4605 
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NOTE—The combined solution is obtained by a combination of all solutions weighted by the probability for the Galactic model (Gal.Mod.) 
and exp( — Ax? /2). 
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cause there is no conclusion about how to combine degen- 
erate solutions. We note that the exp(—Ax?/2) probabil- 
ity might be suffered from systematic errors of the observed 
data. However, the weight from Ay? only has minor effects 
on the “combined results”. Except for KMT-2016-BLG- 
1105 the degenerate solutions have similar physical inter- 
pretations and except for OGLE-2017-BLG-1806 the Ax? 
is small, but for OGLE-2017-BLG-1806 the *combined re- 
sults” are already dominated by the “Close A" solutions due 
to their Galactic-model likelihoods. Due to similar reasons, 
whether to include the phase-space factors also has a minor 
impact on the “combined results". 


5.2. KMT-2017-BLG-1194 


The corresponding CMD shown in Figure 13 is con- 
structed from the OGLE-III field stars within 240" cen- 
tered on the event. The baseline object has (V,J)pase = 
(21.343 + 0.085,19.608 + 0.051), yielding a blend of 
(V — LI)g = (2.15 + 0.39, 20.45 + 0.14). We display 
the blend on the CMD. The source position measured by the 
difference imaging analysis is displaced from the baseline 
object by A0(N, E) = (—26,41) mas. We estimate the 
error of the baseline position by the fractional astrometric 
error being equal to the fractional photometric error (Jung 
et al. 2020), which yields cast = 0.3907 FWHM = 20 
mas. We note that the astrometric error should be underesti- 
mated due to the mottled background from unresolved stars 
and other systematic errors, but the whole astrometric error 
should be not more than twice our estimate. Thus, the base- 
line object is astrometrically consistent with the source and 
the lens within 2c. The blend does not have a useful color 
constraint. We adopt the 3c upper limit of the blended light, 
IL limit = 20.03, as the upper limit of the lens brightness. 

As given in Table 13, the preferred host star is an M dwarf 
located in the Galactic disk, and the planet is probably a 
super-Earth beyond the snow line of the lens system (assum- 
ing a snow line radius asų = 2.7(M/Mo) au, Kennedy & 
Kenyon 2008). 


5.3. KMT-2017-BLG-0428 


The corresponding CMD shown in Figure 13 consists 
of the OGLE-III field stars within 150" centered on the 
event. The baseline object on the CFHT images has Ibase = 
20.056 + 0.063, with an astrometric offset of AO(N, E) = 
(6, —2) mas and an astrometric error of Cass ~ 5 mas. 
Thus, the baseline object is astrometrically consistent with 
the source at about 1c. Because the CFHT images do not 
contain color information, we do not display the blend on the 
CMD. We also adopt the 3c upper limit of the blended light, 
IL limit = 20.81, as the upper limit of the lens brightness. 

As shown in Table 13, the Bayesian analysis indicates an- 
other cold super-Earth orbiting an M dwarf. 


5.4. KMT-2019-BLG-1806 


The CMD of this event is constructed from the OGLE- 
III field stars within 150" centered on the event, shown in 
Figure 13. The baseline object on the KMTC images has 
(V, L)base = (20.155 + 0.125, 18.685 + 0.076). We plot the 
blend on the CMD and find that the blend probably belongs 
to the foreground main-sequence branch and thus could be 
the lens. However, the astrometric offset is A0(N, E) — 
(433, —76) mas and A0(N, E) — (416, —96) mas on the 
CFHT and KMTC images, respectively, so the majority of 
the blended light is unrelated to the lens. We adopt the me- 
dian value of the blended light, Tt limit = 18.8, as the upper 
limit of the lens brightness. 

The results of the Bayesian analysis are given in Table 13. 
The planet is another cold super-Earth, and the preferred host 
is a K dwarf. 


5.5. KMT-2017-BLG-1003 


We use the OGLE-III field stars within 180" centered on 
the event to build the CMD. Combining the measured p from 
the light-curve analysis, we obtain 0g = 0.180 + 0.041 mas 
for the “Outer” solution and 0g > 0.14 mas (3o) for the “In- 
ner" solution. The KMTNet baseline object has (V, I) pase = 
(20.968 + 0.046, 18.780 + 0.028), corresponding to a blend 
of (V — I, I)g = (2.54 + 0.20, 19.83 + 0.10), and we dis- 
play the blend on the CMD. The source-baseline astrometric 
offset is A0(N, E) = (—64, —77) mas, with an astrometric 
error of Cast ~ 12 mas, implying that most of the blend light 
should be unrelated to the event. We adopt the median value 
of the blended light, Jt limit = 19.83, as the upper limit of 
the lens brightness. 

The Bayesian analysis shows that the host star is probably 
an M dwarf located in the Galactic bulge. Again, the pre- 
ferred planet is a cold super-Earth. 


5.6. KMT-2019-BLG-1367 


In Figure 13, we display the position of the source on the 
CMD of stars within 180" around the source. On the KMTC 
images, there is no star within 1."4 around the source posi- 
tion. We thus adopt the detection limit of the KMTC images, 
I — 21.0, as the upper limit of the baseline brightness, yield- 
ing the 3c upper limit of the blended light, JL limit = 21.6. 
Applying Equations (19) and (20) and assuming Dr, < 8 
kpc, this flux constraint corresponds to an upper limit of the 
lens mass of 0.6Mo. 

As shown in Table 13, the Bayesian estimate shows an- 
other cold super-Earth orbiting an M dwarf. 


5.7. OGLE-2017-BLG-1806 


The CMD of this event is constructed from KMTC field 
stars within a 300" square centered on the event position. 
The baseline object, (V, I) base = (22.300 + 0.308, 20.042 + 


0.128), is displaced from the source by 835 mas. Thus, most 
of the blend light should be unrelated to the event. We do not 
show the blend on the CMD and adopt the median value of 
the blended light, Jt, jimit = 20.5, as the upper limit of the 
lens brightness. 

The results of the Bayesian analysis are presented in Ta- 
ble 13, and all solutions indicate a cold super-Earth orbiting 
a low-mass star. The constraints on 7g,, from the light- 
curve analysis are useful. The “Wide” solution has a rela- 
tively large 0g, with a 2c lower limit of 0.60 mas and the 
best-fit value of ~ 1.1 mas, so the corresponding lens system 
is located in the Galactic disk. Then, the “Wide (ug > 0)" 
solution has zg,, < 0 and thus a lens velocity in Galactic 
coordinates of v; ~ 100 km s7}, so this solution is strongly 
disfavored. For the two “Close” solutions, both the mg, < 0 
solutions are slightly disfavored and have relatively higher 
probabilities of a bulge lens system. 

For the “Wide” solution, the predicted apparent magni- 
tude of the lens system is fainter than the source by ~ 0.8 
mag in the H band. In the case of OGLE-2012-BLG-0950, 
the source and the lens have roughly equal brightness and 
were resolved by the Keck AO imaging and the HST imaging 
when they were separated by about 34 mas (Bhattacharya 
et al. 2018). For OGLE-2017-BLG-1806, we estimate that 
resolving the lens and source probably requires a separation 
of 45 mas for the “Wide” solution. We note that the proper 
motions of the two “Close” solutions are ~ 2 masyr- !. If 
high-resolution observations resolve the lens and the source 
and find that pre (e.g., ~ 5 mas yr™ 1) is much higher than 
that of the “Close” solutions, the three solutions can be dis- 
tinguished. Such observations can be taken in 2026 or earlier. 


5.8. KMT-2016-BLG-1105 


To collect enough red-giant stars to determine the centroid 
of the red giant clump, the CMD of this event shown in Fig- 
ure 13 contains KMTC field stars within a 280" x 300" rect- 
angle region. Because the event lies about 80" from the edge 
of the CCD chip, it is displaced from the center of the rect- 
angle region by about 70". The V-band data have insuffi- 
cient signal-to-noise ratio to determine the source color, so 
we adopt the method of Bennett et al. (2008) to estimate the 
source color. We first calibrate the CMD of Holtzman et al. 
(1998) HST observations to the KMTC CMD using the cen- 
troids of red giant clumps. We then estimate the source color 
by taking the color of the HST field stars whose brightness 
are within the 5c of the source star. 

The baseline object has Ibase = 20.729 + 0.125 with- 
out color information, so we do not plot the blend on the 
CMD. The source-baseline astrometric offset is A0(N, E) — 
(73, 166) mas, at about 3c. Because the baseline object is 
marginally detected on the KMTC images, we adopt the me- 
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dian value of the blended light, It limit = 21.7, as the upper 
limit of the lens brightness. 

The Bayesian analysis indicates that the host star is proba- 
bly an M dwarf. Due to a factor of ~ 13 differences within 
the mass ratios of the five degenerate solutions, there is a 
wide range for the planetary mass, from sub-Earth-mass to 
sub-Neptune-mass. Because no solution has a very different 
proper motion from other solutions, future high-resolution 
observations cannot break the degeneracy. However, such 
observations are still important because the measurements 
of the host brightness can yield the host mass and distance, 
which could be used for studying the relation between the 
planetary occurrence rate and the host properties. For the 
“Wide A” and “Close” solutions, the predicted apparent mag- 
nitude of the lens system is fainter than the source by ~ 2 
mag and ~ 3 mag in the H and I bands, respectively. In 
2025, the lens and the source will be separated by = 50 mas 
and may be resolved. 


6. DISCUSSION 


In this paper, we have presented the analysis of seven 
q < 1074 planets. Together with 17 already published and 
three that will be published elsewhere, the KMTNet Anoma- 
lyFinder algorithm has found 27 events that can be fit by 
q < 1074 models from 2016-2019 KMTNet data. For the 
analysis above and in other published papers, all of the local 
minima are investigated, but here for each planet, we only 
consider the models with Ax? < 10 compared to the best-fit 
model. 

Table 14 presents the event name, log q, s, uo, discovery 
method, Ax? compared to the best-fit models, whether it 
has a caustic crossing, anomaly type (bump or dip), and the 
KMTNet fields (prime or sub-prime) of each planet, ranked- 
ordered by logq of the best-fit models. Of them, 15 were 
solely detected using AnomalyFinder, and 12 were first dis- 
covered from by-eye searches and then recovered by Anoma- 
lyFinder, which illustrates the importance of systematic plan- 
etary anomaly searches in finding low mass-ratio microlens- 
ing planets. The seasonal distribution, (5, 8, 8, 6) for 2016— 
2019, is consistent with normal Poisson variations. 

Among the 27 planets, four have alternative possible mod- 
els with q > 1074, and 23 are secure q < 1074 planets. 
Because the detection of q < 1074 planets is one of the ma- 
jor scientific goals of the ongoing KMTNet survey and fu- 
ture space-based microlensing projects (Penny et al. 2019; 
Ge et al. 2022; Yan & Zhu 2022), it is worthwhile to review 
the properties of the 27 planetary events and study how to 
detect more such planets. 


6.1. The Missing Planetary Caustics Problem 


As illustrated by Zang et al. (2021b), the motivation for 
building the KMTNet AnomalyFinder algorithm is to ex- 
hume the buried signatures of “missing planetary caustics” 
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Figure 14. logq vs. log s distribution for the 27 planetary events 
with q < 1074 shown in Table 14, adapted from Figure 11 of Yee 
et al. (2021). The red points represent planets that were solely de- 
tected by AnomalyFinder, and the black points represent planets 
that were first discovered from by-eye searches and then recovered 
by AnomalyFinder. Solutions are considered to be “unique” (filled 
points) if there are no competing solutions within Ax? « 10. Oth- 
erwise, they are shown by open circles. The event KMT-2016-BLG- 
1105 has five degenerate solutions, but we only plot the best-fit 
s > lands < 1 solutions for simplicity. For two solutions that 
are subject to the uo > 0 and uo « 0 degeneracy, we show them 
as one solution and take the average values. The two green dashed 
lines indicate the boundaries for *near-resonant" caustics (Dominik 
1999). 


in the KMTNet data. Zhu et al. (2014) predicted that ~ 50% 
of the KMTNet q < 1074 planets should be detected by 
caustics outside of the near-resonant (Dominik 1999; Yee 
et al. 2021) range. Below we follow the definitions of Zang 
et al. (2021b) and refer to caustics inside and outside of 
the near-resonant range as near-resonant caustics and pure- 
planetary caustics. Contrary to the prediction of Zhu et al. 
(2014), before the application of AnomalyFinder only two 
of ten q < 10 ^^ KMTNet planets were detected by pure- 
planetary caustics. The two cases are OGLE-2017-BLG- 
0173Lb (Hwang et al. 2018a) and KMT-2016-BLG-0212Lb 
(Hwang et al. 2018b). Hence, it is necessary to check the 
caustic types for the planetary sample of AnomalyFinder. 
Figure 14 shows the log q versus log s plot for the 27 plan- 
ets. The red and black points represent planets that were 
first discovered using AnomalyFinder and by-eye searches, 
respectively. The two green dashed lines indicate the bound- 
aries for the near-resonant range. A striking feature is that 
in constrast to the locations of the by-eye planets, of the 


15 AnomalyFinder planets 11 have pure-planetary caustics, 
two have both pure-planetary and near-resonant caustics, and 
only two are fully located inside the near-resonant range. 
In total, at least 13 planets were detected by pure-planetary 
caustics. Thus, the caustic types of the AnomalyFinder plan- 
etary sample agree with the expectation of Zhu et al. (2014), 
and the missing planetary caustics problem has been solved 
by the systematic planetary anomaly search. 


6.2. Caustic Crossing and Anomaly Type 


Zhu et al. (2014) predicted that about half of the KMT- 
Net planets will be detected by caustic-crossing anomalies. 
Jung et al. (2022) found that 16/33 of 2018 KMTNet Anoma- 
lyFinder planets have caustic-crossing anomalies. As shown 
in Table 14, 14/27 of the q < 10-4 planets have caustic- 
crossing anomalies, in good agreement with the expectation 
of Zhu et al. (2014). Thus, the ~ 50% probability of caustic- 
crossing anomalies is likely applicable down to q ~ 1075. 

Zang et al. (2021b) and Hwang et al. (2022) applied the 
AnomalyFinder algorithm to 2018-2019 KMTNet prime- 
field events and found seven newly discovered q < 2 x 1074 
planets. Among them, only OGLE-2019-BLG-1053Lb has 
a bump-type anomaly and the other six planets were de- 
tected by dip-type anomalies. Thus, it is necessary to check 
whether dip-type anomalies dominate the detection of low- 
q planets. As presented in Table 14, the ratio of bump-type 
to dip-type anomalies for the q < 1074 planets is 15 to 12, 
so the two types of anomalies play roughly equal roles in 
the low-q detection. However, of the 12 dip-type anoma- 
lies, nine were solely detected by AnomalyFinder, includ- 
ing eight non-caustic-crossing anomalies. KMT-2018-BLG- 
1988 (Han et al. 2022a) is the only case that the anomaly is 
a non-caustic-crossing dip and was first discovered from by- 
eye searches. Unlike the dip-type anomalies, the four non- 
caustic-crossing bumps were all first noticed from by-eye 
searches. Hence, by-eye searches have proved to be quite in- 
sensitive to non-caustic-crossing dip-type anomalies for low- 
q planets. 


6.3. A Desert of High-magnification Planetary Signals 


Zang et al. (2021b) suggested that the missing plane- 
tary caustics problem was caused by the way that mod- 
elers searched for planetary signatures. Because high- 
magnification events are intrinsically more sensitive to plan- 
ets (Griest & Safizadeh 1998), by-eye searches paid more 
attention to them, while pure-planetary caustics are mainly 
detected in low-magnification events. If this hypothesis is 
correct, we expect that by-eye planets and AnomalyFinder 
planets will have different |uo| and [Uanom| distributions. The 
log |uanom| versus log |uo| distribution of Figure 15 confirms 
our expectation. Except for the two planets that were de- 
tected by pure-planetary caustics, all the other by-eye plan- 
ets, which are located inside the near-resonant range, were 
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Figure 15. log |Wanom| vs. log |uo| distribution for the 27 planetary 
events with q < 1074 shown in Table 14. Colors are the same as the 
colors of Figure 14. Circles and triangles represent prime-field and 
sub-prime-field planets, respectively. The grey dashed line indicates 
[uo| = [tanom]. 


detected with |wo| < 0.05 and |tanom| SS 0.07. The roughly 
one-dex gap of the by-eye planets, at 0.05 < |uo| S; 0.62 
and 0.07 S |uanom| S 0.78, is filled by the AnomalyFinder 
planets?. 

However, there is no planet located at the left lower corner 
of Figure 15, with |uojimit| = 0.0060 and [uanom,limit| = 
0.0158. Although six of the planets were detected in high- 
magnification events (|ug| < 0.01), all the planetary sig- 
nals occurred on the low- and median-magnification regions. 
This desert of high-magnification planetary signals could be 
caused by the insufficient observing cadences of the current 
KMTNet survey. High-magnification planetary signals for 
q < 1074 events are weak and thus require dense observa- 
tions over the peak. There are three known q < 1074 events 
whose planetary signals occurred on the high-magnification 
regions (|uanom| < 0.01). They are OGLE-2005-BLG-169 
with tanom = 0.0012 (Gould et al. 2006), KMT-2021-BLG- 
0171 with [uanom| = 0.0066 (Yang et al. 2022), and KMT- 
2022-BLG-0440 with |Uanom| = 0.0041 (Zhang et al. in 


5 Although it might seem that the correlation could be with anomaly 
brightness rather than |uo| (because smaller |uo| implies a more highly mag- 
nified event), Jung et al. (2022) showed that there is no correlation with event 
brightness at the time of the anomaly between by-eye vs. AnomalyFinder 
detections. On the other hand, Hwang et al. (2022) and Zang et al. (2022a) 
have shown that AnomalyFinder is much better at finding anomalies with 
smaller Ax?. 
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prep). The follow-up data played decisive roles in these de- 
tections and the combined cadences of survey and follow-up 
data are higher than 30 hr 5, while the highest cadence of the 
current KMTNet survey is 8 hr! for about 0.4 deg? from the 
overlap of two T = 4hr7’ fields. 

However, we note that AnomalyFinder used the KMTNet 
end-of-year pipeline light curves, for which the photomet- 
ric quality is not as good as that of TLC re-reductions. For 
the three follow-up planets, the planetary signals only have 
AI < 0.05 mag. Thus, TLC re-reductions may be needed to 
recover such weak signals in the KMTNet data, and we can- 
not rule out the possibility that the desert may also be due to 
the imperfect KMTNet photometric quality. Each year there 
are about 20 events with |uo| < 0.01 observed by KMTNet 
with T > 4 hr !. The current KMTNet quasi-automated 
TLC re-reductions pipeline takes «1 hr of human effort for 
each event (H. Yang et al. in prep), so an optimized sys- 
tematic search for q < 1074 planets in the KMTNet high- 
magnification events can be done very quickly. This search 
could have important implications for future space-based mi- 
crolensing projects, because their tentative cadences are sim- 
ilar to or lower than T = 4 hr! (Penny et al. 2019; Ge 
et al. 2022; Yan & Zhu 2022). If this search demonstrates 
that high-magnification events need denser observations to 
capture the weak planetary signals for low-q planets, one 
could consider conducting (if feasible) ground-based follow- 
up projects for high-magnification events that are discovered 
by space-based telescopes. We also note that for the 2018 
AnomalyFinder planets (Gould et al. 2022; Jung et al. 2022) 
and 2019 prime-field AnomalyFinder planets (Zang et al. 
20222), which are complete now, only one q > 10~* planet, 
KMT-2019-BLG-1953Lb, has |Uanom| < |Uanom,limit|. Fu- 
ture analysis of all the 2016-2019 KMTNet should check 
whether the desert is obvious for more massive planets. 


6.4. Prime and sub-Prime Fields 


In its 2015 commissioning season, KMTNet observed four 
fields at a cadence of T = 6 hr +. To support the 2016-2019 
Spitzer microlensing campaign (Gould et al. 2013, 2014a, 
2015a,b, 2016, 2018) and find more planets, KMTNet mon- 
itored a wider area, with a total of (3, 7, 11, 3) fields at ca- 
dences of T ~ (4, 1, 0.4, 0.2) hr !. The three fields with the 
highest cadence are the KMTNet prime fields and the other 
21 are the KMTNet sub-prime fields. See Figure 12 of Kim 
et al. (2018a) for the field placement. As shown in Table 14 
and Figure 15, the prime fields played the main role in the de- 
tection of q < 1074 planets, as predicted by Henderson et al. 
(2014), and 17 of 27 planets were detected therein. How- 
ever, the sub-prime fields are also important and six of the 
ten lowest-q planets were discovered therein. 

For the six planets with |uo| « 0.01, there is a clear bias in 
cadences, and only one of them was detected from the sub- 


22 ZANG ET AL. 


prime fields. For the prime and sub-prime fields, the cur- 
rent detection rates are 1.25 and 0.25 per year, respectively. 
Because ~ 60% of the KMTNet microlensing events are lo- 
cated in the sub-prime fields, if the sub-prime-field events 
with |uo| < 0.01 can had the same cadence as the prime-field 
events from follow-up observations, each year there would 
be (1.25 x (6096/4096) — 0.25) = 1.6 more q < 107+ 
planets. Because follow-up observations can have higher ca- 
dences and capture the high-magnification planetary signals 
(e.g., Yang et al. 2022), the yield of a follow-up project can 
be at least two q < 1074 planets per year. The reward is not 
only enlarging the low-q planetary sample, but also an inde- 
pendent check of the statistical results from AnomalyFinder 
if the follow-up planets can form a homogeneous statisti- 
cal sample (Gould et al. 2010). However, this reward re- 
quires that the KMTNet alert-finder system should alert new 
events before they reach the high-magnification regions (e.g., 
A > 20). 
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Table 14. Information of 2016-2019 KMTNet AnomalyFinder planetary sample 
with q < 1074 solutions 


Event Name KMTNet Name log q s Juo] Method Ax? Caustic-Crossing Anomaly Type Field 
KB161105! KB161105 —5.194 + 0.248 1.143+0.009 0.171 Discovery 0.0 yes bump sub-prime 
—4.423 + 0.197 1.136 0.011 0.153 2.3 yes 
—4.184 + 0.206 1.106 + 0.013 0.154 2.7 no 
—5.027 + 0.080 0.888 +0.007 0.148 3.5 no 
—4.069 + 0.182 0.892 + 0.005 0.154 4.4 yes 
OB 160007? KB161991 —5.168 + 0.131 2.829 + 0.009 1.253 Discovery yes bump prime 
OB191053? KB191504 —4.885 + 0.035 1.406 40.011 0.373 Discovery yes bump prime 
OB190960* KB191591 —4.830 + 0.041 1.029 0.001 0.0061 Recovery 0.0 yes bump sub-prime 
—4.896 + 0.024 0.997 3 0.001 0.0060 1.0 yes 
—4.896 + 0.024 0.996 + 0.001 0.0059 1.9 yes 
—4.845 + 0.043 1.028 X 0.001 0.0061 2.1 yes 
KB180029° KB180029 —4.737 + 0.047 0.999+0.002 0.027 Recovery 0.0 yes bump sub-prime 
—4.746 + 0.050 1.028+0.002 0.027 0.2 yes 
—4.740 + 0.045 0.999+0.002 0.027 2.1 yes 
—4.736 + 0.050 1.028+0.002 0.027 2.2 yes 
KB191806! KB191806 —4.714 + 0.116 1.035+0.009 0.0255 Discovery 0.0 no dip sub-prime 
—4.717+0.117 1.034 0.009 0.0257 0.4 no 
—4.724 + 0.117 0.938+0.007 0.0260 0.7 no 
—4.734 + 0.109 0.938 + 0.007 0.0251 1.1 no 
OB1701739 KB171707 —4.606 + 0.042 1.540+0.031 0.867 Recovery 0.0 yes bump prime 
—4.195 + 0.068 1.532+0.025 0.844 3.5 yes 
KB171194! KB171194 —4.582 + 0.058 0.806 40.010 0.256 Discovery no dip sub-prime 
KB181988" KB181988 —4.544*0200 0.974 0.03 0.014 Recovery 0.0 no dip sub-prime 
—4.759*0895 — 1.01+0.05 0.014 0.1 no 
KB1908428 KB190842 —4.389 + 0.031 0.983 + 0.013 0.0066 Recovery no bump prime 
KB190253? KB190253 —4.387 + 0.076 1.009+0.009 0.0559 Discovery 0.0 no dip prime 
—4.390 + 0.080 0.929 + 0.007 0.0555 0.3 no 
OB180977? KB180728 —4.382 + 0.045 0.897+0.007 0.147 Discovery yes dip prime 
KB171003! KB171003 —4.373 + 0.144 0.910+0.005 0.179 Discovery 0.0 no dip sub-prime 
—4.260 + 0.152 0.889+0.004 0.179 0.2 no 
OB171806! KB171021 —4.352 + 0.171 0.857 +0.008 0.026 Discovery 0.0 yes bump sub-prime 
—4.392 + 0.180 0.861 + 0.007 0.025 0.2 yes 
—4.441 + 0.168 1.181 0.011 0.026 8.3 yes 
—4.317 + 0.126 1.190 40.012 0.027 8.4 yes 
OB161195!? KB160372 —4.325 + 0.037 0.989+0.004 0.0526 Recovery 0.0 no bump prime 
—4.318 + 0.038 1.079 0.004 0.0526 0.1 no 
OB170448? KB170090 —4.296 + 0.149 3.157 30.022 1.482 Discovery yes bump prime 
—2.705 + 0.045 0.431 +0.004 1.486 5.8 yes 
—3.969 + 0.086 3.593 + 0.045 1.611 9.7 yes 
KB191367! KB191367 —4.303 + 0.118 0.939 + 0.007 0.083 Discovery 0.0 no dip sub-prime 
—4.298 + 0.103 0.976 + 0.007 0.082 0.2 no 
KB170428! KB170428 —4.295 + 0.072 0.882+0.004 0.205 Discovery 0.0 no dip prime 
—4.302 + 0.075 0.915+0.005 0.205 0.1 no 
OB171434!! KB170016 —4.242 + 0.011 00.979 30.001 0.043 Recovery 0.0 yes dip prime 
—4.251+0.012 0.979 30.001 0.043 4.0 yes 
OB181185!? KB181024 —4.163 + 0.014 0.963+0.001 0.0069 Recovery no bump prime 
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OB1811261? KB182064 —4.130 + 0.280 0.852+0.040 0.0083 Discovery 0.0 no dip prime 
—4.260 + 0.290 1.154+0.052 0.0082 2.1 no 

OB180506? KB180835 —4.117 + 0.133 1.059 + 0.021 0.0884 Discovery 0.0 no dip prime 
—4.109 + 0.126 0.861 +0.018 0.0884 0.4 no 

KB181025!4 KB181025 —4.081+0.141 0.937+0.021 0.0071 Recovery 0.0 no bump prime 
—3.789 + 0.133 0.883 + 0.025 0.0086 8.4 no 

OB171691!5 KB170752 —4.013 + 0.152 1.003+0.014 0.0495 Recovery 0.0 yes bump sub-prime 
—4.150 + 0.141 1.058+0.011 0.0483 0.4 yes 

OB180532!9 KB181161 —4.011 + 0.053 1.013 +0.001 0.0082 Recovery 0.0 yes dip prime 
—4.033 + 0.047 1.0112: 0.001 0.0071 2.0 yes 
—3.926 + 0.049 1.013 3-0.001 0.0089 4.6 yes 
—4.016 + 0.076  1.0112-0.001 0.0074 54 yes 

KB160625? KB160625 —3.628 + 0.226 0.741 +0.009 0.073 Discovery 0.0 yes bump prime 
—4.138 + 0.159 1.367 +0.018 0.075 1.0 yes 
—3.746 + 0.291 0.741 +0.009 0.072 1.0 yes 
—4.499 + 0.266 1.358 + 0.015 0.076 3.3 yes 

KB160212*7 KB160212 —1.434 + 0.072 0.829 + 0.007 0.328 Recovery 0.0 yes bump prime 
—4.310 + 0.070 1.427 40.014 0.615 6.6 yes 
—4.315 + 0.099 1.434 + 0.012 0.619 8.0 yes 
—4.082 + 0.080 1.430 + 0.015 0.617 8.7 yes 


NOTE: For each planet, we only consider the models that have Ax? < 10 compared to the best-fit model. “Discovery” represents that the planet was 
discovered using AnomlyFinder, and “Recovery” means that the planet was first discovered from by-eye searches and then recovered by AnomlyFinder. 
Reference: 1. This work; 2. in prep; 3. Zang et al. (2021b); 4. Yee et al. (2021); 5. Gould et al. (2020), Zhang et al. in prep; 

6. Hwang et al. (2018a); 7. Han et al. (2022a); 8. Jung et al. (2020); 9. Hwang et al. (2022); 10. Shvartzvald et al. (2017), Bond et al. (2017), 
Zhang et al. in prep; 11. Udalski et al. (2018); 12. Kondo et al. (2021); 13. Gould et al. (2022); 14. Han et al. (2021); 15. Han et al. (2022b); 

16. Ryu et al. (2020); 17. Hwang et al. (2018b). 


